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Abstract — We consider single-antenna interference networks 
where M sources, each with an average transmit power of P/M, 
communicate with A/ destinations over frequency-selective chan- 
nels (with L taps each) and each destination has perfect knowledge 
of its channels from each of the sources. Assuming that there ex- 
ist error-free non-interfering broadcast feedback links from each 
destination to all the nodes (i.e., sources and destinations) in the 
network, we show that naive interference alignment, in conjunc- 
tion with vector quantization of the impulse response coefficients 
according to the scheme proposed in MukkaviUi et al. , IEEE Trans. 
IT, 2003, achieves full spatial multiplexing gain of M/2, provided 
that the number of feedback bits broadcast by each destination is 
at least M{L - 1) log P. 

I. Introduction 

Cadambe and Jafar [1] proposed a transmission scheme, 
called interference alignment, for single-antenna interference 
networks operating over time-selectiveQ channels and showed 
that this scheme achieves full spatial multiplexing gain. This 
result depends, however, critically on the assumption of each 
source and each destination knowing all the channels in the net- 
work perfectly. In this paper, we show that full spatial multiplex- 
ing gain is achievable even with partial channel state informa- 
tion (CSI) at the sources and the destinations, obtained through 
limited capacity (error-free) broadcast feedback links. In par- 
ticular, we consider an interference network where M single- 
antenna source-destination pairs, denoted by — 
1 , . . . , Af , communicate concurrently and in the same frequency 
band ovex frequency-selective channels with L taps eachrlEach 
source has an average transmit power of P/M and every des- 
tination has perfect knowledge of its channels from each of the 
sources. Our main contribution is to show that naive interfer- 
ence alignment based on vector-quantized impulse responses, 
employing the vector quantization scheme proposed for single- 
user beamforming in [2], achieves full spatial multiplexing 
gain of Af/2, provided that each destination can broadcast at 
the rate M{L — l)logP to all the sources and destinations 
in the network. On a conceptual level, this result shows that 
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'We use the terms time-selective and frequency-selective to denote channels 
that are selective only in time and only in frequency, respectively. 

^Interference alignment, as introduced in [1], does not distinguish between 
time and frequency dimensions. Therefore, although the scheme was originally 
developed for time-selective channels, it can equally well be employed for 
frequency-selective channels. We do not consider time-selective channels as 
vector quantization of the channel coefficients in such channels would require 
non-causal feedback. 



rather than aligning interference perfectly by creating com- 
pletely interference-free signal space dimensions, it suffices to 
ensure that, as the SNR increases, the interference power in 
these dimensions remains bounded. 

Notation: The superscripts ^, ^, and * stand for trans- 
position, Hermitian transpose, and element-wise conjugate, re- 
spectively. CN{Q, cr"^) denotes a circularly symmetric complex 
normal distribution with variance a^. Vectors and matrices are 
set in lower-case and upper-case bold-face letters, respectively. 
I |x| I is the Euclidean norm of the complex vector x and |a;| is the 
absolute value of the complex scalar x. E[-] denotes the expec- 
tation operator. C^^*^ is the set of complex matrices with N 
rows and M columns. The inner product of two column vectors 
a and b of equal dimension is a^b. The diagonal matrix of 
size N X N, with diagonal entries ai, a2, . . . , oat, is denoted by 
diagjoi, 02, ... , a^}. Square brackets [•] and circular brackets 
(•) are used to designate discrete-time and discrete-frequency 
index, respectively. A o B is the Hadamard (or element-wise) 
product of the matrices A and B. log(-) stands for logarithm to 
the base 2 and j = \/— T. The discrete Fourier transform (DFT) 
of the A^-point sequence a[7i], n — 0, . . . , N ~ 1, is defined as 
Tr{a[n]} ^ {1/VN) a[n]e-^-2-w. 

II. System Model 

The L-tap impulse response of the frequency-selective single- 
input single-output (SISO) channel between source Sk,k — 
1 , . . . , M, and destination I?; , i = 1 , • • • , M, is given by 
hi.k[l],l ~ 0, ■ ■ ■ , L~l. The channel coefficients hik[l] remain 
constant throughout the time interval of interest (outage setting) 
and are drawn independently (across i,k,l) from a single con- 
tinuous probability density function such that < |/ii.A;H| < 
oo,yi,k,l, with probability 1. We use a cyclic signal model 
(such as in orthogonal frequency division multiplexing) to con- 
vert the frequency-selective channel Sk — > Vij^k^i, into 
(with N ^ L) parallel frequency- flat channels with coefficient 
hi.k{r), r = 0, . . . , N — 1, for the r-th tone. The input-output 
relation between Si and Vi, for the r-th tone, is then given bjj^ 

y^ir) = ft*,(r)x,(r) + ^ /i* fe(r)a;fc(r) + z,(r) (1) 

k^i 

V ' 

interference 

where yi{r) is the symbol received at destination T>i, Xk{r) de- 
notes the transmit symbol for source Sk and Zi{r) is CAf{0, No) 

^We conjugate the channel coefficients for notational simplicity later on. 



noise at Pj, all for the r-th tone. Defining 

y. = [y.(0) 2/.(i) ... y^{N-l)f 



= [z.(0) z,(l) ... z,{N- 
Hi k = diag{ft.i,fe(0),/ii,fc(l), . 



1)]^ 

.,h,,k{N~l)} 



the input-output relation ([T]i can be rewritten as 



(2) 



E[\xk{r)\' 



< 



fc = l,...,M, r = 0, ...,iV-l. (3) 



The transmit signals obey the power constraints 

P 

M' 

Finally, we shall also need the channel vector ^ and the 
normalized channel vector Wi^k corresponding to the link Sk 
2?„definedasVfc = [/i,,fe[0]'/i,,fc[l] ••• Kk[L-l]f e C^""^ 
andwi^fe = hi,fc/||hi_fe|| e C-^^\ respectively. 

We assume that each destination Vi knows its channels from 
each of the sources Sk perfectly, that is, Vi knows fc,Vfc. 
There exist dedicated non-interfering error- free broadcast feed- 
back links from each destination Vi to all the other terminals in 
the network, that is, to the sources iSfc , Vfc, and to the destinations 
Vk ,ky^i.ln the remainder of the paper, we distinguish between 
a channel feedback phase during which Nf bits of feedback 
are broadcast by each destination and a data transmission phase 
following the channel feedback phase. The channels, being de- 
terministic, are fed back only once during the entire time interval 
of interest so that the transmission rate loss due to the channel 
feedback phase can be assumed to be negligible. Denoting the 
rate of communication for the source-destination pair Si — > "Di 
by Ri, and letting i?sum = 12iLi we say that full spatial 
multiplexing gain is achieved if 

i?sum M 



lim 

P->oo 



logP 



Recall that the spatial multiplexing gain in M source- 
destination pair single-antenna interference networks is upper- 
bounded by M/2 (see [1, Th. 1]). 

III. Interference alignment with perfect CSI at 

ALL NODES 

We next briefly review the concept of interference alignment 
(lA) by adapting the main results of [1] to our setup. 

Each source and each destination knows all the channels in 
the network perfectly. Each source Sk transmits a linear combi- 
nation of dk scalar symbols, x]., x^, . . . , x'l'' , in N frequency 
slots by modulating the symbols onto the transmit direction 



vectors vj., v^, , 



that is. 



'k ! 



m— 1 



•<Nxl 



with I 



,m|l2 



l,andE 



(5) 



,m|21 



where G C,v™ 
P/{Mdk),yk,m. SettingQ = (M-l)(A/-2)-l, the number 



of data symbols dk (corresponding to Sk) and the number of 
tones N are chosen according to (see [1, Appendix III]) 



\tQ, k = 2,3,.. 



(6) 
(7) 



where t is an auxiliary variabl^and the choice of Si to transmit 
{t + symbols in N frequency slots, in contrast to t'^ sym- 
bols for the other sources, is without loss of generality. Each 
destination Vi computes the projections of its received signal 
y,: onto di receive direction vectors u,J , uf , . . . , u^' resulting in 
a total of J2iLi d-i effective input-output relations given by 

(uD^y, = (uD^H^vrrr™ + J2 (O^'H^vf < 



interference 



^^(ur)^H5v^<+(ur)^z. (8) 



interference 



for TO = l,...,d^,i = 1, . . . , M, where uf* e C^''^ with 
||u™|p — 1, Vi, TO. Choosing x™,Vfc, TO, to be i.i.d. Gaussian, 
treating the two interference terms in ([H) as additional noise, 
and assuming that Vi knows the effective channel coefficient 
(u™)-^Hf^v™, Vto, perfectly, the rate of communication over 
the link Si — > Vi is lower-bounded according to 



Md. 



with 



|(u'")^H^-vP| 



k^i p=l 



(9) 
(10) 

(11) 



Each source Sk computes, based on its channel knowledge, 
transmit direction vectors v™, 771 = 1, . . . ,dk, and each desti- 
nation Vi computes, based on its channel knowledge, receive 
direction vectors u™ , m = 1, . . . ,di, that together satisfy the 
following three sets of conditions: 



KuD^Hf^v™! > 0, 

(uD^Hf.vf = 0, 



Vi, m 

Vz, Vm ^ p 
Vfc ^ i,\/m,p 



(12) 
(13) 
(14) 



with the constant c independent of P. It then follows that 1 = 
2^1,2 = 0,Vi, and the spatial multiplexing gain achieved by lA 



'*We employ the auxiliary variable t, partly to simplify our exposition, and 
partly to keep our presentation consistent with [1]. The precise role of t will 
become clear later. 



is lower-bounded according to 

M di / 
r, E E log 1 . 

iim > iim 

P^oo log P P^oc N log P 



' 2 



iV 



that is, full spatial multiplexing gain, in the sense of Q, is 
achieved. Under the assumption of every node (i.e., every source 
and every destination) knowing all the channels in the network 
perfectly, one way to find vectors u™ , satisfying ([T2jl-([T4l) 
is provided in [1, Appendix III]. The basic idea is that each 
Sk computes, based on its knowledge of all the channels in 
the network, a set of linearly independent transmit direction 
vectors , v| , . . . , v)^*" such that all the vectors corresponding 
to interference from Sk,k ^ i, at Vi (that is, the vectors 
H^^v^, Vfc i,p = 1, . . . , dfc) span an {N — di)-dimensional 
complex subspace of C^. Consequently, d, dimensions remain 
completely interference-free. Each Vi, in turn, computes, based 
on its knowledge of all the channels in the network, a set of di 
unit-norm receive direction vectors , u| , . . . , u^' that spans 
the d^-dimensional interference-free subspace corresponding to 
the link Si Vi, thereby satisfying ( [T4] i. Moreover, it was 
shown in [1, Appendix III] that if the vectors fif^jV™, Vm, along 
with the vectors H^j.v^,Vfc ^ j,Vp, span C^, then Vi can 
choose the di receive direction vectors u™, Vm, such that along 
with ( fT4l ), both ([T2| and ( [T3] l are satisfied as well. It turns out 
that, in the frequency-selective case, this is possible provided 
that L > {{t+ l)Q - l)/{3tQ) (the proof of this statement is 
similar to [3, Th. 6.4] and the details are provided in [4]). 

The developments in the remainder of this paper are based 
on the simple observation that if the interference power terms 
1; 1 and Ti,2, for all i, are not equal to zero, but upper-bounded 
by a constanj^iWe/^enc/enf of P, full spatial multiplexing gain 
is still achieved. The key to realizing this will be a vector 
quantization scheme, which satisfies ( [T2] ) and ensures that both 
KuD^H^vf |2,Vz,Vm ^ and KuD^H^v^p, Vfc ^ 
Vm,p, scale as 1/P when P ^ oo. It will turn out that the 
vector quantization scheme developed in [2] and [5] for beam- 
forming in single-user frequency-flat multi-input multi-output 
(MIMO) channels satisfies this condition. 

IV. Interference alignment with limited feedback 

We start by recalling that each destination Vi knows the chan- 
nel coefficient vectors h^jt, Vfc, (and hence, Wi.fe, Vfc) peifectly. 
Knowledge of 'Wi^k,'^k, at all the sources and all the other 
destinations is obtained through feedback. Specifically, each 
Vi broadcasts, during the channel feedback phase, quantized 
versions of Wi.fe, Vfc, to 5fe, Vfc, and Vk,yk ^ i. We shall next 
describe the vector quantization and feedback scheme used. 



^To be precise, full spatial multiplexing gain is achieved even if X^.i and 

0. 



Xi,2 scale as a function of P, say f{P), such that limp_,tx> '°iog jT'* 



Relegating the details to [4], we note, however, that this does not result in a 
reduction of the required feedback rate scaling (in P). 



A. The vector quantization and feedback scheme 

The vector quantization scheme works on unit norm vec- 
tors and quantizes the vector Wi.fc e C^^^ to the unit norm 
vector Wi e C^^^ using Nj^ bits. The corresponding quan- 
tizer codebook A therefore contains 2^"* vectors, that is, A = 
{wi, W2, . . . , WjNrf }. The quantization policy is as follows: 



Wi,fc = arg max {|wf Wj^fcl} • 



(15) 



The quantizer codebook is chosen as the solution of the follow- 
ing Grassmannian line-packing problem [6 ]: Find the maximum 
number of unit-magnitude vectors in such that the absolute 
value of the inner product between any two of the vectors is less 
than cos((5), where 5 G (0, 7r/2]. Here, 6 is an auxiliary variable 
whose role will become clear later. 

The line-packing problem: If A^pack is the solution to the line- 
packing problem and {pi, p2, . . . , PATp^^^l is a set of vectors 
corresponding to this solution, we choose 

2^" = A^pack and set 

{wi,W2, ■ • • ,W2iVrf} = {pi,P2,. . . ,PArp^,^J. (16) 

This approach was used in [2] and [5] for beamforming in 
single-user MIMO channels. 

Quantization error: We define the quantization error 

Ad(wi^fc,Wi_fc) as 



Ad(wi^fe,Wi,fc) = y 1 - |w5,.Wi,fep. 
The maximum quantization error is then given by 



= max 
xec-^ Jlxii 



(17) 



where Wx G ^ is the unit-magnitude quantized version of x G 
obtained according to ([TSj. 

We will need an upper bound on A™^^ in terms of A^^. 
While such a bound is known [2], we will provide a derivation, 
partly for completeness, and partly to get the bound in a form 
required for our proof. We start by noting the following two key 
properties of the chosen set of quantization vectors: 

i) The following relation holds between A^^^ and sin((5) (see 



[5, Th. 3]): 



Na < ^ sin((5) 



-2(L-1) 



sin((5) < 2 



92(-L-l) 



(18) 



ii) The maximum quantization eiTor A™^'' is upper-bounded 
by sin((5) . This can be proved by contradiction as follows: 

• Let us assume that A™'^'*^ > sin((5). 

• Further, let Xq be the unit-magnitude vector, quan- 
tized to Wx„, that corresponds to the maximum quan- 
tization error, that is. 



Xo = arg max 

xec^ J|x||=i 



1 - Ix^w 



(19) 



Then, we have 



1 - |x?w 



> sin((5) 



< 1 - si-a^iS) = 

< cos((5) 

xf pi| < cos((5), I = 1, 



(20) 

cos^id) (21) 
(22) 

...,2^''. (23) 



• However, ( [23) implies that in the line-packing 
problem, there exists a solution set of vectors 
{Pi: P2, • ■ • , P2"d : Xo}, such that the absolute value 
of the inner product between any two vectors is less 
than cos{d). In other words, 2^"* (and consequently 
Npack) is not the maximum possible number of unit- 
magnitude vectors with the absolute value of the in- 
ner product between any two vectors being less than 
cos((5) and can hence not be the solution of the line- 
packing problem, which results in a contradiction. 
Thus, our premise A™^^ > sin((5) must be incorrect 
and we must necessarily have 

Ar'' < sm{6). (24) 
Inserting ( |24] l into ( [T8] l, we get the desired upper bound: 



1 



(25) 

\2-2(L-l) ) 

Number of feedback bits: During the channel feedback phase, 
each destination Vi broadcasts Nd bits for the realization 
Wi^fe, k = 1, . . . , M, to all sources and destinations (except to 
itself, of course) in the network, resulting in a total of 



Nf = MNd bits 



(26) 



being broadcast by each destination. Each source Sk therefore 
receives a total of M'^N^ bits of (error-free) feedback from 
all the destinations and each destination Vi receives a total of 
(Af — \)MNd bits from the destinations I?fe,fc ^ i. Each 
source and each destination can therefore recreate the quantized 
normalized channel vectors Wi_fc,Vi, k, and the key to proving 
the main result of this paper is to determine a value of Nd that 
ensures the achievability of full spatial multiplexing gain with 
naive lA based on Wj^^:, Vi, k. 

B. Transmission scheme and achievability of full spatial multi- 
plexing gain 

Each source and each destination first converts its received L- 
dimensional quantized vectors \v'i^i,,yi, k, into A^-dimensional 
vectors [m^k[0] •■• Wi,kW - 1]]^ through zero- 

padding. It then computes the TV-point DFTs 

m,k{r)^J'r{m,kM}^ r = 0,...,iV-l (27) 
and organizes the results into the quantized channel matrices 

W,,fc = diag{u-,,fc(0),m,,fe(l), . . .,w,AN - 1)}. (28) 



lA is now performed naively assuming that H 



i.k 



w 



i,k : 



Vi, k, that is, each source Sk computes its transmit direction 
vectors v™, m = 1, . . . , dk, and each destination Vi computes 



its receive direction vectors u™ , m 



. ,di, from 



Wijt,Vi,fc (rather than from Hi^k)- Sk transmits a linear 
combination of dk scalar symbols, x^, . . . , a;^*" , in N 
frequency-slots by modulating the symbols onto the vectors 



Xfc 



(29) 



771 — 1 

where e C, v™ e 



:^><iwith||v^"i|2 



l,andE[|5 



.771121 



P/ {Mdk), Vfc, m. The number of data symbols dk (correspond- 
ing to Sk) and the number of tones N are chosen according 
to (|6]) and (|7|, respectively. Each destination Vi computes the 
projections of its received signal onto the receive direction 
vectors u^^ , u^, . . . , u^' resulting in a total of X]f=i '^i effective 
input-output relations given by 



Y^Y.^K)"^?,kK< + io"'^ 

k^i p —1 



Vi, m (30) 



-•JVxl 



with ||u™||2 = l,Vi,m. Defining h, k = 

/^,,.(iV-l)Fandb:"/^(u^)*ov^, 



where u™ e 

[/iz.fc(O) Kk{l) ... h,k{N-l)V andbT;:^ ^ 
we can rewrite the input-output relations ( [jo] i as 

(uD^y, = hf,b: 



;_ff ^777,771^771 



+ EE^^5bD'< + ("r)''^- Vz,m. (31) 

k^i p—1 

Choosing x^yi^m, to be i.i.d. Gaussian, treating the two 
interference terms in ( (3T| as additional noise and assuming 
that the destination Vi knows the effective channel coefficients 
hf'^b™'™,Vm, perfectly, the rate of communication over the 
link Si 



Vi is then lower-bounded according to 

Mdi \ 



with 




(32) 



k^ip —1 

Recall that in lA with perfect CSI, the conditions (fT2|)-(fl4l) 
are satisfied. Defining = (u 

are equivalent to 



o v^, these conditions 



Ih^-b"'™ 



I > c > 0, 



Vi, m 



(33) 



h^^f = 0,Vi,Vm ^ p, and h^^kKk'' = 0,Vfc 7^ i,Vm,p, 
respectively. Naive lA entails finding vectors u™ , satisfying 
the following conditions: 

> c> 0, yi,m (34) 

\li,\/m 7^ p (35) 



1,1 1,1 



^i,k"i.k 



0, 



(36) 



where w^^fe = [wi^k{0) Wi,fc(l) • • ■ Wj,fc(^^ - 1)]^- As noted 
earlier, one way to find vectors u™ , satisfying ([34|-([36| is 
provided in [ 1 , Appendix III] . The key point here is that although 
based on imperfect CSI, this choice of u™, vf , Vi, m, p, results 
in full spatial multiplexing gain and, in addition, this can be 
realized with a feedback rate of M{L — 1) log P. We proceed 
with the proof of this statement. 

Since |lwi^i|| — ||wi.i|| ~ 1 (from Parseval's theorem) 
and Wi.ijb™'*' e C^^^, V«, to,p, we can always find vectors 
qi,q2, . . . ,q7v-2 such that {w,,i, b™'P/||b"'''|| , qi, q2, . . . , 
qAf-2},?Ti ^ p, is an orthonormal basis for C^. Expanding 
hj i into this orthonormal basis, we get 



|h,,,;||2 = |hf,W,; 



N-2 



E MM 

1=1 

Vi, Vm ^ p 



which yields 
P 



Mdi 



h^b™-'' 



< 



P 



Mdi 



< 



P 



< 



Mdi 
4P 
Jldi 



ib:7ip (iih. 



p 

Mdi 


|bI7in|h.,,f 


(^^ 




1 


P 








1 


Mdi 


l|h.,,;|| 





bi7iH|h,,f (Ar'^)2 




(37) 
(38) 

(39) 
(40) 

yi^m^p (41) 



where (|39| follows from Parseval's theorem and ( |4T] ) is obtained 
by invoking ( |25] l. 

A similar expansion of fc into the orthonormal basis 
{w,,fc,b"//||b"'/||,qi,q2, . . . , qAr-2}, A: 7^ i, yields 



Mdk 



2 4P - ,,- ,,9 




Vfc ^ i,\fm,p. (42) 
If we now choose Nd ~ (L — 1) log P, we get from ( |4T] i that 



P 

mI, 



2,2 Z.Z 



2 4P - 

< ||b'".P||2 Jj. 



4||br-nP ||h. 



1 



" V ' 

=l/p 

Vi, Vm 7^ p 



and from ((42| that 

P 2 

"■i,k"i,k 



Mdk 




which implies that, with Nd = {L — l)logP, the overall 
interference power in the rate lower bound ( (32] i is upper- 
bounded by a constant independent of P according to 



,k ■ 



(43) 



p^ra k^i p — 1 

Finally, we note that as P ^ oo , the vector quantizer codebook 
size 2^'' also tends to infinity and the maximum quantization 
error tends to zero according to ( p5] l. The resolution of the 
vector quantizer therefore becomes arbitrarily high and we 
obtain Wi^i \ii^J\\\ii,i\\,wh\ch implies Wi_j hi^i/||hi_j||. 
Substituting this into ([34]| yields a condition equivalent to ( [33] l, 
that is, 

h^-b"'" 



Ih^b-'^l 



> c > 0, 
Vi, m. 



Vi, TO 



(44) 



>||h,,,||c>0 Vi,TO. (45) 

Consequently, using ( |43| ), the spatial multiplexing gain achieved 
by naive lA is lower-bounded according to 

Psiim 



lim 

P^cc log P 



log 1 



M di 

^EE 

i^l rn — 1 

N 



lim 

P^oo 



iVlogP 
(from ^ and (j?!) 



(46) 



which proves that full spatial multiplexing gain is achieved. We 
complete the proof by noting that the number of bits fed back 
(broadcast) by each destination for achievability of full spatial 
multiplexing gain using naive lA is given, according to ( |26] l, by 

Nf = MNd = M{L - 1) log P. 
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